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I.  Overview 

Graphene  is  an  attractive  material  for  spintronics  due  to  its  high  mobility  and  the  low  intrinsic 
spin-orbit  and  hyperfme  coupling,  which  should  lead  to  excellent  spin  transport  properties.  In 
2007,  graphene  became  the  first  material  to  exhibit  gate  tunable  spin  transport  and  spin 
precession  at  room  temperature  [12],  However,  the  spin  injection  efficiency  was  low  and  the  spin 
lifetime  was  much  shorter  than  predicted  theoretically.  These  issues  defined  the  immediate 
challenges  for  the  field.  The  original  goals  of  this  project  were  to  investigate  two  broad 
directions  in  graphene  spintronics:  (1)  Spin  injection  and  transport  in  graphene:  significantly 
improve  the  characteristics  of  spin  injection  and  spin  lifetime;  explore  novel  geometries  such  as 
pn  junctions  and  nanoribbons;  investigate  spin  transport  in  epitaxial  graphene.  (2)  Spin  transport 
in  doped  graphene:  investigate  the  effects  of  surface  impurity  scattering  on  charge  and  spin 
transport;  find  evidence  for  magnetism  in  doped  graphene  via  exchange  fields;  generate  spin- 
Hall  effect  generated  by  spin-orbit  fields.  We  have  made  substantial  progress  in  both  of  these 
areas  including  some  breakthrough  results  that  has  established  our  group  as  among  the  top  two 
groups  in  graphene  spintronics. 

(1)  Spin  Injection  and  Transport  in  Graphene 

We  have  made  important  progress  in  spin  injection  and  transport  [1,2,3,7,9,10]  including: 
Tunnelins  spin  injection  into  2ravhene.  The  low  spin  injection  efficiency  into  graphene  is  due  to 
the  conductivity  mismatch  between  the  ferromagnetic  metal  (Co)  and  the  single  layer  graphene 
(SLG).  To  alleviate  this  problem  and  enhance  the  spin  injection  efficiency,  we  grew  atomically 
smooth  MgO  tunnel  barrier  with  Ti02  seed  layer.  With  tunneling  contacts,  the  non-local  spin 
signal  was  found  to  be  as  high  as  130  Q  at  room  temperature,  with  a  spin  injection  efficiency  of 
30%.  This  is  currently  the  highest  spin  injection  efficiency  observed  in  graphene  spin  valves  [7]. 
Lon2  spin  lifetimes  in  srayhene.  Enhanced  spin  lifetimes  are  observed  using  tunneling  contacts 
[7].  This  indicates  that  the  short  spin  lifetimes  reported  before  are  due  to  the  contact- induced  spin 
relaxation  from  the  ferromagnetic  electrodes.  Using  tunneling  contacts,  we  observed  the  spin 
lifetimes  as  high  at  771  ps  in  single-layer  graphene  (SLG),  1.0  ns  at  4  K  in  SLG,  and  6.2  ns  at  20 
K  in  bilayer  graphene  (BLG)  [9].  These  are  currently  the  highest  values  reported  in  the  literature. 
At  low  temperatures,  contrasting  behaviors  of  gate  voltage  dependence  of  the  spin  lifetime  are 
observed  between  SLG  and  BLG,  which  suggest  different  mechanisms  for  spin  relaxation  in 
SLG  and  BLG  [9]. 

(2)  Spin  and  Charge  Transport  in  Doped  Graphene 

Due  to  the  extreme  surface  sensitivity  of  graphene,  there  is  a  unique  opportunity  to  modify  the 
spin-dependent  properties  of  graphene  through  a  proximity  effect.  By  introducing  adatoms  or 
thin  films  on  top  of  graphene,  one  could  induce  exchange  fields,  tunable  spin-orbit  coupling,  or 
hyperfme  fields.  These  will  be  extremely  useful  for  manipulation  of  spin  in  graphene.  Our 
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studies  utilize  in  situ  transport  measurements,  which  allow  adsorbates  to  be  deposited  in 
ultrahigh  vacuum  and  transport  measurements  to  be  taken  without  exposing  the  sample  to  air. 

We  have  several  results  [4, 5, 6, 8]  including: 

Charsed  impurity  scatterins  and  char2e  transfer  from  transition  metal  impurities:  We 
performed  the  first  in  situ  study  of  transition  metal  doping  on  graphene.  For  adsorbates  deposited 
onto  graphene  at  low  temperature  to  keep  the  adatoms  separated  (i.e.  to  avoid  aggregation),  we 
found  that  all  materials  studied  including  Fe,  Ti,  Au,  and  Pt  donated  electrons  to  graphene.  This 
was  especially  surprising  for  Pt,  which  has  a  bulk  workfunction  of  ~6.0  eV  that  is  much  larger 
than  graphene  (~4.5  eV).  In  such  a  case,  one  would  expect  the  Pt  to  act  as  an  electron  acceptor. 
This  shows  for  atomic  adsorbates  (or  small  clusters),  the  simple  consideration  of  bulk 
workfunction  is  insufficient  to  not  applicable  [4].  We  also  investigated  the  role  of  clustering  [5] 
and  the  role  of  insulating  vs  metallic  character  of  the  adsorbates  [8]. 

Role  of  charsed  impurity  scatterins  on  spin  relaxation:  By  investigating  the  effect  of  gold 
dopants  on  spin  transport,  we  concluded  that  charged  impurity  scattering  is  not  the  dominant 
source  of  spin  relaxation  in  graphene,  (ref  6)  This  is  a  very  important  result  because  the  issue  of 
spin  relaxation  is  among  the  most  important  in  graphene  spintronics  due  to  the  huge  discrepancy 
between  the  long  spin  lifetimes  expected  theoretically  (microseconds)  and  shorter  lifetimes 
observed  experimentally  (hundred  picoseconds).  Several  possible  mechanisms  have  been 
proposed  for  the  spin  relaxation.  This  work  eliminates  charged  impurity  scattering  as  one  of  the 
possibilities. 


II.  Detailed  Description  of  Results 


(1)  Spin  Injection  and  Transport  in  Graphene 


Graphene 


Fig.l  :  Schematic  drawing  of  a  graphene  spin  valve  in  the  non-local  spin  transport 
geometry. 


All  experiments  in  this  section  were  based  on  nonlocal  spin  transport  in  graphene  spin  valves 
consisting  of  an  ultrathin  sheet  of  graphene  (single  or  bilayer)  contacted  by  ferromagnetic  cobalt 
electrodes  (Figure  1).  The  devices  were  connected  in  the  non-local  geometry,  with  a  current 
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source  to  inject  spins  into  the  graphene  and  a  voltmeter  to  measure  the  spin  diffusion  away  from 
the  current  loop.  A  pure  spin  current  (without  net  charge  current)  flows  from  injector  to  detector. 
The  spin  transport  was  measured  in  the  non-local  magnetoresistance  geometry,  where  a  magnetic 
field  was  applied  in-plane  along  the  long  axes  of  the  Co  electrodes.  Spin  lifetimes  were  measured 
in  the  Hanle  geometry,  where  a  magnetic  field  was  applied  out-of-plane  to  generate  spin 
precession  of  the  pure  spin  current.  The  Si  substrate  was  used  as  a  backgate  to  change  the  carrier 
concentration  of  the  graphene. 

At  the  start  of  the  grant  period,  the  two  overriding  issues  in  graphene  spintronics  were  the 
short  spin  lifetimes  and  the  poor  spin  injection  efficiency.  One  route  to  increasing  the  spin 
injection  efficiency  was  to  utilize  a  tunnel  barrier  to  alleviate  the  conductance  mismatch  problem 
for  spin  injection.  Building  on  our  earlier  results  of  growing  atomically  flat  MgO  on  graphene 
via  submonolayer  Ti02  seed  layer,  we  became  the  first  to  demonstrate  tunneling  spin  injection 
into  graphene. 

Figure  2  shows  the  electrical  properties  of  graphene  spin  valves  with  Co  electrodes  and 
Co/MgO/Ti02/graphene  tunnel  junctions.  Figure  2a  is  the  bias  dependence  of  the  differential 
contact  resistance  (dV/dl  of  the  contact).  The  peak  near  zero  bias  indicates  the  non-linearity  of 
the  I-V  curve  of  the  junction,  which  is  characteristic  of  tunneling.  Figure  2b  shows  the 
temperature  dependence  of  the  zero  bias  differential  contact  resistance.  The  slight  decrease  as  a 
function  of  increasing  temperature  also  indicates  that  the  transport  across  the  junction  is  via 
tunneling  as  opposed  to  pinholes.  Together,  these  characteristics  show  that  transport  across  the 
junctions  is  by  tunneling. 


Fig.  2  :  Electrical  characterization  of  contacts,  (a)  Bias  dependence  of  differential  contact 
resistance  shows  that  the  I-V  curve  is  nonlinear,  (b)  Temperature  dependence  of  zero-bias 
contact  resistance. 


Investigating  the  spin  injection  and  transport  in  graphene  spin  valves  with  tunneling  contacts 
proved  to  be  very  interesting  [7].  Figure  3  compares  the  non-local  spin  signal  for  devices  with 
transparent  contacts  (Figure  3a)  and  tunneling  contacts  (Figure  3b).  Interestingly,  for  devices 
with  similar  distance  between  spin  injection  and  spin  detector  (~2  microns),  the  non-local  spin 
signal  was  2,000  larger  for  the  tunneling  contact  than  the  transparent  contact.  The  spin  injection 
efficiency,  which  is  the  spin  polarization  of  the  injected  carriers,  was  ~30%  for  tunneling 
contacts  while  only  ~1%  for  the  transparent  contacts.  This  manifested  the  expected  enhancement 
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of  spin  injection  efficiency  by  using  tunnel  barriers  to  alleviate  the  conductance  mismatch 
problem. 


(a)  T  ransparent  Contacts 


(b)  Tunneling  Contacts 
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Fig.  3  :  Enhancement  of  spin  injection  efficiency  by  tunneling  contacts,  (a)  Non-local  spin 
transport  measurement  for  graphene  spin  valve  with  transparent  contacts,  (b)  Non-local 
spin  transport  measurement  for  graphene  spin  valve  with  transparent  contacts. 


Investigating  the  effect  of  tunnel  barriers  on  spin  lifetime  provided  a  surprising  result  [7]. 
Figure  4  compares  the  spin  lifetime  measure  by  Hanle  measurement  for  transparent,  pinhole,  and 
tunneling  contacts.  The  enhancement  of  spin  lifetime  up  to  771  ps  revealed  contact-induced 
effects  as  the  origin  of  the  strongest  spin  relaxation  in  graphene  spin  valves.  The  contact  effects 
could  be  interfacial  spin-flip  scattering,  dephasing  by  inhomogeneous  fringe  fields,  or  escape 
time  effects. 


(a)  Transparent  contacts  (b)  Pinhole  contacts  (c)  Tunneling  contacts 
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Fig.  4  :  Enhancement  of  spin  lifetime  due  to  tunneling  contacts,  (a-c)  Hanle  measurement 
of  spin  lifetime  in  graphene  spin  valve  with  transparent  contacts  (panel  a),  pinhole 
contacts  (panel  b),  and  tunneling  contacts  (panel  c).  Red  (black)  curve  is  for  parallel 
(antiparallel)  magnetization  alignment  of  injector  and  detector  contacts. 
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Fig.  5:  Relationship  between  spin  lifetime  and  diffusion  coefficient  in  graphene,  (a)  As  a 
function  of  gate  voltage  (carrier  concentration),  single  layer  graphene  exhibits  a  similar 
behavior  for  spin  lifetime  and  diffusion  coefficient,  which  is  characteristic  of  Elliot- Yafet 
spin  relaxation,  (b)  As  a  function  of  gate  voltage  (carrier  concentration),  bilayer  graphene 
exhibits  opposite  behavior  for  spin  lifetime  and  diffusion  coefficient,  which  is 
characteristic  of  Dyakonov-Perel  spin  relaxation. 


Using  high  quality  tunneling  contacts,  we  investigated  spin  relaxation  in  single  layer 
graphene  (SLG)  and  bilayer  graphene  (BLG)  [9].  At  low  temperatures,  we  found  that  the  spin 
lifetimes  in  BLG  were  typically  longer  than  SLG.  The  longest  reported  value  for  BLG  was  6.2 
ns.  Also  of  strong  interest  was  comparing  the  relationship  between  spin  relaxation  and 
momentum  scattering  (via  diffusion  coefficient  D).  In  SLG,  the  spin  lifetime  Ts  and  diffusion 
coefficient  D  were  found  to  have  a  similar  dependence  on  carrier  concentration  (Figure  5a).  This 
suggests  an  Elliot- Yafet-like  mechanism  for  spin  relaxation  (spin  relaxation  generated  by 
momentum  scattering).  On  the  other  hand,  in  BLG  the  spin  lifetime  Xg  and  diffusion  coefficient  D 
were  found  to  have  opposite  dependence  on  carrier  concentration  (Figure  5b).  This  suggests  a 
Dyakonov-Perel-like  mechanism  for  spin  relaxation  (spin  relaxation  generated  by  precession 
from  internal  spin-orbit  fields,  which  is  suppressed  by  momentum  scattering).  Some  recent 
theoretical  studies  based  on  Rashba  spin-orbit  domains  provide  a  possible  explanation  for  the 
coexistence  of  Elliot- Yafet-like  and  Dyakonov-Perel-like  properties  [13]. 

We  investigated  the  role  of  bias  dependence  on  spin  injection  and  spin  extraction  in  SLG 
spin  valves  with  transparent  contacts  [3].  The  bias  dependence  is  a  topic  of  high  interest  within 
spintronics  for  semiconductor  and  metal  systems.  The  unique  aspect  of  graphene  is  that  the 
polarity  of  the  carriers  can  be  tuned  from  electrons  to  holes.  Specifically,  the  graphene  band 
structure  exhibits  an  electron-hole  symmetry  that  does  not  exist  in  conventional  semiconductors. 
In  our  experiments,  we  observed  an  asymmetry  in  the  spin  injection/extraction  for  electrons  and 
holes.  When  the  carriers  were  holes  and  the  bias  was  for  spin  extraction,  the  non-local  spin  signal 
was  reduced  (Figure  6).  This  did  not  occur  for  spin  extraction  by  electrons,  implying  a  breaking 
of  the  electron-hole  symmetry.  We  believe  this  broken  symmetry  is  due  to  the  n-type  doping  of 
graphene  by  the  direct  Co  contact. 
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Fig.  6  :  Bias  dependence  of  non-local  spin  signal  for  electrons  (solid  dots)  and  hole  (open 
dots).  For  these  two  samples,  the  spin  injection  efficiency  exhibits  a  sharp  drop  for  holes 
under  negative  bias.  This  corresponds  to  hole  spin  “extraction”,  where  the  holes  are 
driven  from  the  graphene  into  the  Co  electrode. 
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(2)  Spin  and  Charge  Transport  in  Doped  Graphene 

We  developed  an  ultrahigh  vaeuum  ehamber  that  allows  for  the  deeoration  of  graphene 
surfaees  using  moleeular  beam  epitaxy  (MBE)  effusion  eells,  in  situ  variable  temperature 
magnetotransport  measurement  (15  -  300  K,  0-1000  Oersted  field),  and  sample  transfer.  This 
enabled  studies  of  adsorbed  atoms  on  the  surfaee  of  graphene  and  their  effeet  on  eharge  and  spin 
transport,  all  without  exposing  the  sample  to  air.  The  in  situ  approaeh  was  first  used  by  the 
Fuhrer  group  to  study  the  effeets  of  potassium  doping  on  eharge  transport  in  graphene  [14]. 
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Fig.  7:  Effect  of  transition  metal  doping  |  < 

on  graphene  charge  transport,  (a)  SEM 
image  of  device,  (b)  AFM  image  of 
transition  metal  clusters,  (c-e)  Evolution  a 
gate-dependent  conductivity  as  Ti,  Fe, 
and  Pt  atoms  are  deposited  onto  graphene. 

(f)  Dirac  point  shift  as  a  function  of  coverage  for  Ti, 
Fe,  and  Pt.  (g-i)  Mobility  change  for  electrons  and 
holes  for  Ti,  Fe,  and  Pt,  respectively. 
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Our  initial  investigation  centered  on  the  effects  of  transition  metal  doping  on  charge  transport 
in  graphene  [4].  To  our  knowledge,  this  was  the  first  study  of  this  type  using  transition  metals. 
Figure  7a-7c  shows  the  changes  in  gate-dependent  conductivity  as  Ti,  Fe,  or  Pt  atoms  were 
deposited  onto  the  graphene  surface  at  room  temperature.  In  all  cases,  the  Dirac  point  shifted  to 
the  left,  indicating  that  the  transition  metals  were  donating  electrons  to  the  graphene.  This  was 
especially  surprising  for  Pt,  which  has  a  bulk  workfunction  of  ~6.0  eV  that  is  much  larger  than 
graphene  (~4.5  eV).  In  such  a  case,  one  would  expect  the  Pt  to  act  as  an  electron  acceptor.  This 
showed  for  atomic  adsorbates  or  small  clusters,  the  simple  consideration  of  bulk  workfunction 
cannot  explain  the  observed  behavior.  The  second  effect  was  that  the  mobility  was  reduced  by 
the  transition  metal  dopants,  as  indicated  by  reduction  of  slope  of  the  conductivity  vs.  gate 
curves.  A  summary  of  the  shift  in  Dirac  point  and  reduction  of  mobility  are  shown  in  Figures  7f- 
7i.  Figure  7f  shows  results  for  several  samples.  While  there  were  sample-to-sample  variations, 
the  main  trend  was  that  Ti  produced  the  strongest  electronic  doping,  while  Pt  produced  the 
weakest  n-type  doping.  This  is  consistent  with  the  trend  of  the  bulk  workfunction.  Figures  7g-7i 
show  the  decrease  in  mobility  for  the  electrons  and  holes  for  Ti,  Fe,  and  Pt,  respectively.  The  Ti 
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Mobility,  /y  (10^  cmWs) 


exhibited  the  fastest  reduction  in  mobility,  which  is  consistent  with  the  strongest  doping  (Fig.  If) 
assuming  charged  impurity  scattering.  Likewise,  Pt  showed  the  slowest  reduction  in  mobility, 
which  is  consistent  with  the  weakest  doping  (Fig.  If).  Interestingly  Fe  showed  an  asymmetry  of 
the  electron  and  hole  mobilities,  which  might  be  due  to  hybridization  of  the  Fe  with  graphene. 
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Fig.  8  :  Effect  of  cluster  formation  on  doping  and  scattering  by  gold  adatoms,  (a)  Gate- 
dependent  conductivity  as  gold  atoms  are  deposited  onto  the  graphene  surface  at  low 
temperature  (18  K).  (b)  Mobility  at  low  temperature  as  a  function  of  gold  coverage,  (c) 
Relation  between  Dirac  point  shift  and  inverse  mobility,  showing  a  power  law  of  1.3 
which  is  consistent  with  pointlike  charged  impurity  scattering,  (d-f)  The  Dirac  point, 
mobility,  and  temperature,  respectively,  as  a  function  of  time  as  the  sample  is  heated  up. 
The  final  data  points  (open  circles)  are  the  values  after  a  second  cooldown,  which  shows 
that  the  changes  are  irreversible. 


We  investigated  the  effect  of  clustering  of  adsorbates  by  utilizing  low  temperature  deposition 
of  gold  adatoms  [5],  which  are  known  to  have  a  large  surface  mobility  at  room  temperature.  The 
low  temperature  deposition  freezes  the  positions  of  the  gold  atoms  and  subsequent  annealing  to 
room  temperature  allows  the  atoms  to  move  and  form  clusters.  Figure  8a  shows  the  initial 
deposition  at  low  temperature,  which  exhibited  strong  electronic  doping  (Dirac  point  shift)  with 
nearly  one  electron  per  gold  atom,  similar  to  behavior  seen  for  low  temperature  potassium 
doping  [14].  Analyzing  the  relationship  between  Dirac  point  shift  and  inverse  mobility  showed  a 
power  law  consistent  with  point-like  Coulomb  scatterers  (Figure  8c).  Figure  8d-8f  shows  the 
Dirac  point  shift  and  mobility  as  the  sample  was  warmed  up  to  room  temperature.  With 
increasing  cluster  size,  the  Dirac  point  shift  was  reduced  and  the  mobility  increased.  This 
behavior  is  consistent  with  theoretical  predictions  of  Katsnelson,  Guinea,  and  Geim  [15]. 
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Fig.  9:  Comparing  metallic  vs.  insulating  adsorbates  on  graphene,  (a-c)  Evolution  of 
gate-dependence  conductivity  during  metallic  Ti  deposition  (panel  a),  cluster  formation 
(panel  b),  and  oxidation  of  Ti  in  oxygen  gas  to  form  insulating  TiOi  (panel  c).  (d) 
Evolution  of  mobility  during  this  process,  (e)  Evolution  of  conductivity  minimum  during 
this  process. 


We  also  compared  the  effects  of  metallic  vs.  insulating  adsorbates  on  charge  transport  in 
graphene  [8].  The  strategy  was  to  deposit  metallic  titanium  adsorbates  and  then  post-oxide  to 
convert  the  adsorbates  to  insulating  Ti02.  In  this  manner,  the  comparison  between  metallic  vs. 
insulating  was  done  for  the  same  spatial  distribution  of  adsorbates  on  the  same  piece  of 
graphene.  Figure  9  shows  the  gate-dependent  conductivity  with  the  intial  deposition  of  Ti  atoms 
(Figure  9a),  cluster  formation  (Figure  9b),  and  subsequent  post-oxidation  (Figure  9c).  The  entire 
process  was  performed  at  room  temperature.  Figure  9d  shows  the  evolution  of  the  mobility, 
whereas  the  Figure  9e  shows  the  evolution  of  the  conductivity  minima  during  this  process.  Both 
quantities  are  plotted  against  the  Dirac  point  shift.  Converting  metallic  Ti  adsorbates  to 
insulating  Ti02  by  post-oxidation  generally  reversed  the  effect  of  the  adsorbates.  As  the  Ti  was 
oxidized,  the  mobility  increased  as  the  Dirac  point  came  back  toward  its  initial  value. 
Interestingly,  the  mobility  extrapolated  to  a  value  lower  than  the  initial  mobility  (Fig.  9d).  On  the 
other  hand,  the  conductivity  minimum  extrapolated  up  to  its  initial  value  (Fig.  9e).  Because 
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conducitity  minima  depends  on  the  electron-hole  puddles  what  arise  from  Coulomb  interaction 
with  charged  impurities,  while  the  mobility  can  be  affected  by  both  the  charged  impurity  and 
short  range  scattering.  The  result  of  these  extrapolations  was  to  show  that  the  TiOi  adsorbates 
indeed  contribute  to  short  range  scattering. 
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Fig.  10:  Effect  of  Au  adatom  deposition  on  spin  lifetime  and  diffusion  coefficient  (inset) 
for  electrons  (red),  holes  (blue),  and  at  the  Dirac  point  (black). 

We  also  utilized  the  in  situ  approach  to  investigate  the  mechanism  of  spin  relaxation  in 
graphene  [6],  one  of  the  major  topics  in  graphene  spintronics.  To  investigate  the  role  of  charged 
impurity  scattering  on  spin  relaxation,  we  utilized  gold  adsorbates  deposited  at  low  temperature. 
From  earlier  studies  (Figure  8),  low  temperature  gold  was  found  to  act  as  point-like  charged 
impurity  scatterers.  In  addition,  the  gold-graphene  chemical  interaction  was  expected  to  be  weak. 
The  fact  that  gold  could  introduce  additional  spin-orbit  coupling  was  not  ideal,  but  gold 
remained  as  the  best  choice  for  charge  impurity  scattering  of  all  materials  that  we  had  examined 
up  to  that  point.  Figure  10  shows  the  main  result,  which  was  that  the  addition  of  charged 
impurity  scattering  from  gold  did  not  reduce  the  spin  lifetime,  even  though  substantially  reduced 
the  mobility.  In  fact,  the  spin  lifetime  exhibited  a  slight  increase.  The  lack  of  a  major  decrease  in 
spin  lifetime  proved  that  spins  are  robust  against  charged  impurity  scattering  [6]. 
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